Microscopic Theory for the Quantum to Classical Crossover in Chaotic Transport 
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We present a semiclassical theory for the scattering matrix S of a chaotic ballistic cavity at finite 
Ehrenfest time. Using a phase-space representation coupled with a multi-bounce expansion, we 
show how the Liouville conservation of phase-space volume decomposes S as S = S cl © 5 qm . The 
short-time, classical contribution <S cl generates deterministic transmission eigenvalues T = or 1, 
while quantum ergodicity is recovered within the subspace corresponding to the long-time, stochastic 
contribution 5 qm . This provides a microscopic foundation for the two-phase fluid model, in which 
the cavity acts like a classical and a quantum cavity in parallel, and explains recent numerical data 
showing the breakdown of universality in quantum chaotic transport in the deep semiclassical limit. 
We show that the Fano factor of the shot-noise power vanishes in this limit, while weak localization 
remains universal. 
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PACS numbers: 73.23.-b, 74.40.+k, 05.45.Mt 

Quantum mechanics strongly alters the classical the- 
ory of transport. In mcsoscopic systems, spectacular new 
phenomena of purely quantum origin emerge, such as 
weak localization, universal conductance fluctuations Q , 
and sub-Poissonian shot-noise Q- In both mesoscopic 
disordered diffusive metals and clean chaotic ballistic 
systems, these phenomena have been found to exhibit 
a universality well captured by Random Matrix Theory 
(RMT) 0. For instance, in a time-reversal-symmetric 
chaotic cavity, RMT correctly predicts that weak local- 
ization reduces the conductance below its classical value 
by <5c/ rm t = —1/4, that the variance of the conductance 
ofmtG?) = 1/8 (conductances are expressed in units of 
e 2 /h), and that the shot-noise power is reduced by the 
Fano factor F rmt = 1/4, below the Poissonian value of 
2e(J), (I) being the average current through the sample. 

Despite the many successes of RMT it has become 
clear, following the work of Aleincr and Larkin Q, that 
universality is broken for those quantum chaotic systems 
in which the Ehrenfest time becomes a relevant time 
scale. The Ehrenfest time, te, corresponds to the time 
it takes for a minimal wavepacket to be stretched over 
a length comparable to a classical characteristic length 
of the system. For a closed chaotic cavity of linear size 
L, the Ehrenfest time is t£ ~ A _1 ln[L/Ap], in term of 
the classical cavity's Lyapunov exponent A and the Fermi 
wavelength Af |5| ■ For the same cavity attached to two- 
leads of width W there is a second shorter Ehrenfest time 
t|! = Tg — 2A _1 ln[L/W] 0]. These two quantum mechan- 
ical time scales differ only by a classical quantity. For 
times shorter than the Ehrenfest time, the quantum me- 
chanical time evolution of a narrow wavepacket is given 
by the solution of classical equations of motion [f| 0, 
thus strong deviations from the RMT of transport emerge 
once is no longer negligible. It is for example well es- 
tablished that, as Tg' c — > oo, the Fano factor disappears 
P. 1^. ITol ITlj and parametric conductance fluctuations re- 
main universal while sample-to-sample conductance fluc- 
tuations do no t fill TT3 | . However the analytical predic- 
lEfthal 



tion of Rcfs. weak localization corrections also 



disappear has recently been challenged numerically |l4j . 

These observations are all consistent with the phe- 
nomenological two-phase fluid model [Til Il5j . In this 
model, introduced in Ref. 0, it is assumed that elec- 
trons with short dwell time < are transmitted deter- 
ministically through the cavity, while quantum mechan- 
ical stochasticity is recovered for those with longer dwell 
time > Tg. Accordingly, the system splits into two cavi- 
ties in parallel; a classical, deterministic cavity, with an 
effective number of modes N c \ = N(l — exp[— t^/td]), 
and a quantum mechanical, stochastic cavity having 
N qm = A^cxpf— rg/ro] modes (td is the average dwell 
time). With the further assumption that the stochas- 
tic cavity obeys RMT [U Il5j . the observed numerical 
behavior of shot-noise and of conductance fluctuations 
follows 0, 0, The universal value of paramet- 

ric conductance fluctuations is maintained as long as 
N qm oc h c g +1 ^ XTD 3> 1, where the effective Planck's 
constant h e s = ^f/L (Xf is the Fermi wavelength). 
For fully developed chaotic dynamics Atd S> 1, and so 
this condition for N qm is always satisfied in the semi- 
classical limit of h c e — * ^l|- However the prediction 
that weak localization remains universal (independently 
of ?i e ff) is in direct contradiction of the analytical pre- 
diction Sg oc exp[— t^/td] of Rcfs. 0, 0|. Both this 
controversy and the successes of the phenomenological 
two-phase fluid model call for a microscopic foundation 
of this model. The purpose of this letter is to provide 
such a foundation. 

We consider an open, two-dimensional chaotic quan- 
tum dot, ideally connected to two external leads, each 
carrying N ^> 1 modes. We require that the size 
of the openings to the leads is much smaller than the 
perimeter of the cavity but is still scmiclassically large, 
1 -C N <C L/Xf- This ensures that the chaotic dynamics 
inside the dot has enough time to develop, i.e. Atd 3> 1. 
The classical dynamics within the cavity can be described 
by its Birkhoff map, which provides a recursive relation 
for the transverse momentum p T and position q T of in- 
cidence at the r th collision with the cavity's boundary. 



Quantizing this map leads to a unitary, M x M Floquct 
operator V where M = Jnt[h~^]. A particle initially in a 
state ipo occupies a state ip T = V T V>o after r collisions at 
the cavity's boundary. The transport properties of this 
system derive from its scattering matrix |l6| 

which we write in terms of N x JV transmission (t) and 
reflection (r) matrices. From S, the system's conduc- 
tance is given by g = Tr(t^t). To construct S from V, we 
couple the cavity to the left (L) and right (R) leads by in- 
troducing a 2N x M projection matrix P = P^ + p( R \ 
Expressed in the basis of channel modes, the projection 
matrices read Pnm = 1 if m = n 6 { m i L } and 

Prim = otherwise. The sets {mf 1 ' R '} are the N com- 
ponents of V ideally connected to the modes of the L or 
R lead respectively. The energy-dependent S becomes 

m 

oo 

S(e) = ^exp[i(T + l)e]5 T , (2a) 

S T (e) = P[V(1 - P T P)] T VP T . (2b) 

The T th term in the above Taylor expansion corresponds 
to the time-evolution operator for a particle colliding ex- 
actly t times at the boundary of the cavity before exiting. 

To maximally resolve the underlying classical dynam- 
ics, we wish to express S in an appropriate orthonormal 
basis. We start from a complete set of coherent states 
centered on a discrete von Neumann lattice [l^. While 
this set is not orthogonal, it is nevertheless complete (and 
not ovcrcomplctc) . It is then possible to recursively con- 
struct an orthonormal Phase Space (PS) basis from it, 
{\(p, q)i)h,B.}, where each PS basis state is exponentially 
localized in both position and momentum and covers a 
phase-space area 0(K) [2(j. The algorithm for construct- 
ing this basis is similar to those used in wavelet analy- 
sis , and we will discuss it elsewhere |2l|. Here we 
only show such a PS state in Fig. ^ The transformation 
from lead modes to this orthonormal basis is a unitary 
transformation S — > W SU and hence leaves the trans- 
port properties of iS unaffectedx. We will now show that 
the multi-bounce expansion of Eqs. @, expressed in this 
orthonormal PS basis elegantly connects quantum trans- 
port to the underlying classical dynamics. 

Two cases have to be considered when evaluating the 
matrix elements ((j>, q)i\S T \(p, q)j) of <S T . If \{p,q)j) sits 
on classical trajectories which all leave the cavity at the 
same time and through the same lead, then its quantum 
mechanical time-evolution is given by the classical one. 
This requires first that the classical trajectory starting at 
(p, q)j leaves the cavity after a time r < |{| (this en- 
sures that S T \(p, q)j) is still a well localized wavepackct), 
and second that at the time of escape the wavepackct 
is not located too close to the edge of the lead, and so 
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FIG. 1: Squared amplitude of a PS state (solid line) versus 
the normalized coordinate q/qo, qo — L\/2-KTi c ff_. For compar- 
ison we also show a coherent state (dashed line). 

avoids partial reflection. When these two conditions are 
fulfilled, the whole wavepacket escapes the cavity at a 
single time r through either the L or R lead (not both) 
in a state S T \ (p, q)j } . In all other instances this is not 
the case. For wavepackets that exit the cavity either at 
times t > Tg or after one or more partial reflections at 
the edge of a lead, diffraction and quantum coherence 
play an important role. Thus we split S into a short- 
and a long-time contribution |22| . 

(MilSlMj) = ((p^^ + SniP^h), (3a) 

s cl = J2 s r ; s qm = J2 s t- ( 3b ) 

r<T° T>T° 

The classical Liouville conservation of PS volume guar- 
antees that <S cl and S qm act on different, nonoverlap- 
ping subsets of the PS basis. To see this, we recall that 
the evolution of wavepackets that are initially localized 
in both position and momentum is very well approxi- 
mated semiclassically as a sum over classical trajecto- 
ries 0,E1|. Thus if a classical trajectory contributes to 
«S C , then it does not contribute to S qm . From this Li- 
ouville blocking, we conclude that the PS states divide 
into two non-overlapping subsets (up to exponentially 
small corrections) (i) the fully classical ones which leave 
the cavity after short times r < and are transmit- 
ted or reflected deterministically, and (ii) the quantum 
ones which stay inside the cavity for longer times t > 
and/or are partially reflected and transmitted at collision 
at the edge of a lead |24j • From the dwell time distribu- 
tion p(r) = TpJ 1 cxp[— t/td] in chaotic cavities, one gets a 
number N c \ = N(l — exp[— Te/td]) of states in subset (i) 

and N qm = iVexp[— t^/td] oc ft cff 1+1 ^ ATD 3> 1 of states in 
subset (ii). These considerations lead to our first result 
that S has the block-diagonal decomposition 

S = 5 cl ©5 qm , (4) 

where <5> cl ' qm both have the structure given in Eq. QJ. 
Thus the transmission matrix T = t T t is also block- 
diagonal 

T = T cl ©T qm , (5) 
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with iV cl , qm by AT cl , qm matrices T cl < qm = (t c Um)t t d,qin 
Accordingly, the system's conductance is given by the 
sum of two contributions 



9cl + 5qm, 



(6) 



where g cX = N d /2 and g qm ~ N qm /2. Note that g qm gives 
the contribution to the total conductance carried by the 
stochastic scattering states. As we will see below, it con- 
tains the weak localization corrections. Eqs. 1)416(1 give 
a microscopic justification of the two-phase fluid model 
previously postulated in Refs.[Tll Il5j . where the finite 
t e ' transport properties of a quantum chaotic system 
are given by those of two separated cavities put in par- 
allel. The decomposition (0J of 5 is naturally connected 
to the underlying classical dynamics. In fact, because of 
the fmiteness of td, classical trajectories injected into a 
cavity are naturally grouped into PS transmission and 
reflection bands [13, despite the ergodicity of the asso- 
ciated closed cavity. These bands are best visualized by 
considering PS cross-sections of the left and right leads, 
and plotting the initial and final coordinates of trajecto- 
ries transmitted from left to right. This is sketched in 
Fig. [3] All trajectories within one band follow neighbor- 
ing paths and exit at the same time r through the same 
lead. Because of the chaotic classical dynamics, bands 
with longer escape times r are narrower, having a width 
(and hence a PS area) scaling like oc exp[— At]. The 
Ehrenfest time is the time at which this area becomes 
smaller than h e g, consequently for times longer tha n rg 
no band is big enough for a PS state to fit within it |l5j |. 
The meaning of Eq. Q is now obvious: The matrix iS cl 
acts in the subspace of PS states which are entirely inside 
a single band, while iS qm acts in the subspace of PS states 
which overlap many bands, exiting at different times, and 
in general through different leads. 

We now discuss the properties of T cl and T qm . 
The classical part of the transmission matrix is 
y^/^o fc£t T '. However each PS state in t cl follows a 
band that exits at a well defined time, thus only diago- 
nal terms in the double sum over t, t' are non-zero, and 
one has 



(7) 



We stress that this is not the usual diagonal approxima- 
tion and so we do not require any energy or ensemble 
averaging here. Once T cl is written as in Eq. J7J, it is 
easy to show that it is diagonal. To see this one recalls 
first that for r < rg, S T \(p,q)j) = \{p, q)i) is still a nar- 
row wavepacket, which however has been stretched and 
rotated by the classical dynamics || (this is sketched 
in Fig. I2J). This can be included for each band indi- 
vidually into a r-dependent area-preserving coordinate 
transformation, so that a basis made of these stretched 
PS states is complete and orthonormal within each band 
(relationships between bands can be ignored because Li- 
ouville blocking ensures that bands never overlap). This 




L lead at time 
(injection) 



R lead at time X 
(exit) 



FIG. 2: (color online) Phase space of the injection (Left) 
and exit (Right) leads. The shaded area depicts a classical 
transmission band which goes from L to R in time r. Super- 
imposed on the L lead is a square von Neumann lattice with 
spacing (2Tvh s fi) 1 ^ 2 ■ A PS state is placed at every vertex on 

the lattice, forming a complete orthonormal basis. Each PS 

1/2 

state initially covers a circle of radius 0[h c g ] in phase space. 
For t < Tg, this state is still localized, but is stretched under 
the dynamics, as indicated by the dark ellipse on the R lead. 



means that, after a r-dependent unitary transformation 
of the bra-states, ((p, q)i\t T \(p, q)j) = unless the state 
\(p, q)j) in the L lead is the one state that transmits to 
the (stretched and rotated) state |(p, q)i) in the R lead, 
then ((p,q)i\t T \(p,q)j) = cxp[iifij]. Eq. JJJ shows that 
T cl is unaffected by this transformation, hence 



E 



[trU = Si] XT, i e {l.-.-iVd}, (8) 



with transmission eigenvalues Tj which are either zero or 
one. While this result was anticipated [25| , to the best of 
our knowledge it is derived here for the first time directly 
from a microscopic theory. Thus T cl cannot contribute to 
noise, moreover, all quantum phases (such as tpij) cancel 
in T cl and hence it carries no quantum coherence. Both 
shot-noise and coherent effects like weak localization are 
carried solely by T qm which we now focus on. 

Unlike in T cl , the off-diagonal terms r ^ r' strongly 
influence T qm = Y] tlt T r. These terms give rise to 
coherence and to the pseudo-randomness of transmission 
phases, as many initial PS states may now be partially 
transferred to the same final PS state. Conversely, each 
initial PS state is partially transferred to many final PS 
states 0- Taking the Liouville blocking into account, 
t qm is ergodic within its own subspace (being decou- 
pled from t cl ). Because it is made up of all the long- 
time contributions with t > Tg, it contains dynamically 
diffractive contributions necessary for weak localization 
(Richter-Sieber pairs H||). We further note that its size 



N. 



-1+1/Atd 



qm 



"eft 



cxd in the semiclassical limit. 



This is all we need to generalize the semiclassical the- 
ory of weak localization to finite r^' /to- The approach 
of Ref. [2(| can be followed provided Eq. Q is taken 
into account by considering an effective cavity with a 
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number of modes N qm , a restricted phase-space area 
^qm.(E) = exp[— Tg/r D ] and a dwell time distri- 

bution p qm (r) = 6(r - rg)cxp[-(r - rg)/r D ]/T D . A 
careful analysis of the dwell time for the Richter-Sieber 
pairs (and use of an optimal basis |20| ) yields 5g = 
-1/4(1 + OKXtu)- 1 }) independent of r°' c jl^. Thus 
for good chaotic systems, Atd ^> 1, weak localization 
remains universal in the classical limit. 

Turning our attention to shot-noise, we readily see that 
the reduction of the fraction of quantum channels re- 
produces the predicted exponential suppression of the 
Fano factor in the large rg regime @, F = ^2 n T n (l — 
T n )/T,n T n « N qm /N = exp[-rg/T D ]. We expect there 
will soon be a semiclassical theory to explain the uni- 
versal value F = 1/4. Indeed we believe that a theory 
of transmission along the lines of Ref. [2^ is possible, 
which would show that the properties of S qm are cap- 
tured by one of the circular ensembles of random ma- 
trices. Assuming that this is the case, we have the full 
distribution P(T) of transmission eigenvalues P(T) = 
aP rmt (T) + (l-a) [S(T) + 6(1 - T)] /2, where P rmt (T) = 



7r-i[T(l - TV^}/ 2 and a = exp[-rg/r D ], as postu- 
lated in Ref. [Rj and numerically observed in Ref. [Tl| . If 
confirmed analytically, this would give the full counting 
statistics of the current for any value of t^/td, including 
the correct Fano factor F = (1/4) exp[-rg/r D ] HUll- 

In summary, we have constructed a theory for trans- 
port in quantum chaotic ballistic systems in the regime of 
finite c / t d- Our theory confirms the separation of the 
system into two subsystems, and thus provides a micro- 
scopic foundation for the two-phase fluid model [Til Hb| . 
Weak localization is predicted to remain universal in low- 
est order in 1/ Atd, for any h e s, in agreement with numer- 
ical data 01 > but in contradiction with earlier theories 
0, 0| . Whether S qm is a random matrix remains an 
open question. 
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